Are Galaxy Clusters Suggesting an Accelerating Universe Independent of SNe la 

and Gravity Metric Theory? 
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A kinematic method to access cosmic acceleration based exclusively on the Sunyaev-Zel'dovich 
effect (SZE) and X-ray surface brightness data from galaxy clusters is proposed. By using the 
SZE/X-ray data from 38 galaxy clusters in the redshift range 0.14 < z < 0.89 [Bonament et al., 
Astrop. J. 647, 25 (2006)], we find that the present Universe is accelerating and that the transition 
from an earlier decelerating to a late time accelerating regime occurred relatively recent. Such 
results are fully independent on the validity of any metric gravity theory, the possible matter- 
energy contents filling the Universe, as well as on the SNe type la Hubble diagram from which 
the present acceleration was inferred. The ability of the ongoing Planck satellite mission to obtain 
tighter constraints on the expansion history through SZE/X-ray angular diameters is also discussed. 
Two simple simulations of future Planck data suggest that such technique will be competitive with 
supernova data besides being complementary to it. 
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Introduction. The dimming of distant type la super- 
novae observed by two different group of astronomers one 
decade ago lead to unexpected and landmark conclusion: 
the universal expansion is speeding up and not slowing 
down as beheved since the early days of observational 
cosmology [H 0. 

Such phenomenon is normally interpreted as a dynamic 
influence of some sort of dark energy whose main effect is 
to change the sign of the decelerating parameter g(z)[3'|. 
Another possibility is that the cosmic acceleration is a 
manifestation of new gravitational physics (rather than 
dark energy) that involves a modification of the left hand 
side (geometric sector) of the Einstein field equations. In 
this sort of theory the Friedmann equation is modified 
and a late time accelerating stage is obtained even for 
a Universe filled only with cold dark matter (CDM)[3]. 
At present, the space parameter associated with the cos- 
mic expansion is too degenerate, and, as such, it is not 
possible to decode which mechanism or dark energy com- 
ponent is operating in the cosmic dynamics 

Currently, SNe type la are not only the powerful stan- 
dard candles available but still provides a unique direct 
access to the late time accelerating stage of the Uni- 
verse. Naturally, this a rather uncomfortable situation 
from the observational and theoretical viewpoints even 
considering that ten years later, the main observational 
concerns about errors in SNe type la measurements, like 
host galaxy extinction, intrinsic evolution, possible selec- 
tion bias in the low redshift sample seem to be under 
controlis]. 

A promising estimator fully independent of SNe type 
la and other calibrators of the cosmic distance ladder is 



the angular diameter distance {Da{z)) from a given set of 
distant objects. It has also been recognized that the com- 
bination of SZE [6*1 and X-ray surface brightness measure- 
ments may provide useful angular diameters from galaxy 
clusters 0, i, i, d [Hi . 

On the other hand, since the mechanism causing the 
acceleration is still unknown, it is interesting to inves- 
tigate the potentialities of SZE/X-ray technique from a 
more general viewpoint, that is, independent of the grav- 
ity theory and the matter-energy contents filling the Uni- 
verse. The better strategy available so far is to consider 
the same kind of kinematic approach which has been 
successfully applied for determining the transition decel- 
eration/acceleration in the past by using SNe type la 
measurements [H [H Q [11 [H . 

In this letter, we employ a purely kinematic descrip- 
tion of the universal expansion based on angular diameter 
distances of clusters for two different expansions of the 
deceleration parameter. As we shall see, by using the 
Bonamente et al.^V^ sample we find that a kinematic 
analysis based uniquely on cluster data suggests that the 
Universe undergone a "dynamic phase transition" (de- 
celeration/acceleration) in a redshift z « 0.3. Further, 
it is also shown that the Planck satellite mission data 
must provide very restrictive limits on the space param- 
eter, thereby opening an alternative route for accessing 
the expansion history of the Universe. 

Angular Diameter and Kinematic Approach. Let us 
now assume that the Universe is spatially flat as moti- 
vated by inflation and WMAP measurements [l^ ■ In this 
case, the angular diameter distance in the FRW metric 
is defined by (in our units c = 1), 
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FIG. 1: a) Contours in the qo — qi plane for 38 galaxy clusters data [ij considering q{z) = qo + qiz. The best fit to the pair 
of free parameters is (go, gi) = ( — 1.35, 4.2). For comparison we have shown the straight lines denoting the transitions redshifts 
for two different ACDM models: Zt = 0.9 for JIa ~ 0.8 and zt — 0.15 for JIa ~ 0.43. b) Likelihood function for the transition 
redshift. The best fit is zt = 0.32. 



where H(z) = d/a is the Hubble parameter, and, q{z), 
the deceleration parameter, is defined by 



q{z) EE - — 



dH-\z) 



dt 



1. 



(2) 



In the framework of a flat FRW metric, Eq. (1) is 
an exact expression for the angular diameter distance. 
As one may check, in the case of a linear two-parameter 
expansion for q{z) — qo + zqi, the above integral can 
analytically be represented as 



Da{z) = ^^t/^ ' e^^gi^»-^H7(gi-go,(^ + l)gi) 

-Ho 

-7(91-90,91)], (3) 

where go and 9i are the values of q{z) and its redshift 
derivative, dq/dz evaluated at z = while 7 is an incom- 
plete gamma (18| . By using the above expressions we may 
get information about qo, qi and, therefore, about the 
global behavior of q{z). In principle, a dynamic "phase 
transition" (from decelerating to accelerating) happens 
at q{zt) = 0, or equivalently, zt = —qo/qi. Another inter- 
esting parametrization is q(z) = go + 9i^/(l + z) [l5l.[l^. 
It has the advantage to be well behaved at high redshift 
while the linear approach diverges at the distant past. 
Now, the integral (1) assumes the form: 

Da{z) = ^^^e^^gr^'°+'^^[7(9i +90,91) 
-no 

-7(91 +90, 91/(1 + ^))], 



(4) 



where qi now is the parameter yielding the total cor- 
rection in the distant past {z ^ 0, g(z) = 90 + 91) 
and 7 is again the incomplete gamma function. Note 
that in this case the transition redshift is defined by 
Zt = -90/(90 + 9i)- 



Constraints from Galaxy Clusters. The SZE is a small 
distortion on the Cosmic Microwave Background (CMB) 
spectrum provoked by the inverse Compton scattering of 
the CMB photons passing through a population of hot 
electrons. Observing the temperature decrement of the 
CMB spectrum towards galaxy clusters together the X- 
rays observations, it is possible to break the degeneracy 
between concentration and temperature thereby obtain- 
ing Da{z). Therefore, such distances are fully indepen- 
dent of the one given by the luminosity distance, Dl{z). 

Let us now consider the 38 measurements of angu- 
lar diameter distances from galaxy clusters as obtained 
through SZE/X-ray method by Bonamente and cowork- 
ers In our analysis we use a maximum likelihood 
determined by a statistics 



X {z\v) = 2^ — -2 — : 



Vao.,? 



(5) 



where "DaoA is the observational angular diameter dis- 
tance, (j-Dao i is the uncertainty in the individual distance, 
(^stat is the contribution of the statistical errors added in 
quadrature (« 20%) and the complete set of parameters 
is given by p = (i?o, 9o, 9i)- For the sake consistency, the 
Hubble parameter Hq has been fixed by its best fit value 
= 80km/s/Mpc. 
Linear Parameterization: q — qo + qiz. In Figs. 1(a) 
and 1(b) we show, respectively, the contour in the plane 
qo - qi (68.3%, 90% and 95.4% c.l.) and likelihood of 
the transition redshift from the Bonamente et al. sam- 
ple. The confidence region (la) are —2.6 < qo < —0.25 
and 13 < qi < —3. Such results favor a Universe with 
recent acceleration (qo < 0) and a previous decelerat- 
ing stage [dq/dz > 0). From Fig. 1(a) we see that he 
best fits to the free parameters are qo ~ — 1.35,gi = 4.2 
while for the transition redshift is zt — 0.32 (see Fig. 
1(b)). Note the presence of a forbidden region forming a 
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FIG. 2: a) Contours in the 50 — Qi plane for 38 galaxy clusters data Tlf considering q{z) = qo + qiz/{l + z). The best fit to 
the pair of free parameters is (go, 9i) = (—1.43, 6.18). b) Likelihood function for the transition redshift and the associated best 
fit at zt = 0.30. Comparing with Figs. 1(a) and 1(b) we see that the results are weekly dependent on the parameterizations. 



trapezium. The horizontal hne at the top is defined by 
qi = 0, which leads to an infinite (positive or negative) 
transition redshift. Note also that the segment at 45% 
defines the infinite future (zj = —1). In addition, one 
may conclude that the vertical segment on the left clos- 
ing the trapezium is also unphysical since it is associated 
Zt < —1, thereby demonstrating that the hatched trapez- 
ium is actually a physically forbidden region (for a sim- 
ilar analysis involving luminosity distance see [IBl)- For 
comparison we have also indicated in Fig. 1(a) the tran- 
sition redshifts zt = 0.15 corresponding to a flat ACDM 
with r^A = 0.43, as well as, zt — 0.9 corresponding to 
~ 0.8. 

2nd Parameterization: q = qo + qiz/1 + z. In Figures 
2(a) and 2(b) we display the corresponding plots for the 
second parameterization. The confidence region (la) is 
now defined by: —2.4 < qo < —0.5 and 13.5 < qi < 0. 
Such results also favor a Universe with recent acceleration 
{qo < 0) and a previous decelerating stage {dq/dz > 0). 
As indicated in Fig. 2(a), the best fits to the free pa- 
rameters are go = —1.43 and qi — 6.18 while for the 
transition redshift is a little smaller zt — 0.3 (see Fig. 
1(b)). It should be noticed the presence of the forbidden 
region (trapezium) with a minor difference in comparison 
with Fig. 1(a), namely, as an effect of the parameteri- 
zation, the horizontal line now is at the bottom. Note 
also that a decelerating Universe today (go > 0) is only 
marginally compatible at 2cr of statistical confidence. 

The results in the qo — qi planes for both cases suggest 
that: (i) the Universe had an earlier decelerating stage 
{qi = dq/dz > 0), and (ii) the Universe has been ac- 
celerating (go < 0) since z ~ 0.3. A similar result has 
been previously obtained using SNe type la as standard 
candles by Shapiro and Turner p^. 

Prospects for Planck Satellite Mission. Let us discuss 
the potentiality of the SZE/X-ray technique when future 
data from Planck satellite mission become available [l9|. 



The mission is a project from European Space Agency 
whose frequency channels were carefully chosen for mea- 
surements of thermal Sunyaev-Zeldovich effect. In prin- 
ciple, the Planck satellite will see (through SZE) about 
30,000 galaxy clusters over the whole sky with a signif- 
icant fraction of clusters near or beyond redshift unity. 
However, since accurate angular diameter measurements 
require long SZE/X-ray integrations nobody expects that 
all observed clusters might have useful distance measure- 
ments to constrain cosmological parameters. Therefore, 
it is interesting to simulate two realistic samples of angu- 
lar diameter distances (ADD) by using a fiducial model 
to -Dao"? = DAizi,qo,ql,H*), where H*, q* and ql are 
the best fit values to the linear case obtained from Bona- 
mente et al. sa 

mpleini- 

Table 1 

z range Clusters bins Clusters/bin ADDi Error 

2^0) (P,0) (P,0) 

[0.0,0.5] 100,500 10 10,50 15%, 10% 

[0.5, 1.0] 70, 350 10 7, 35 17%, 12% 

[1.0, 1.5] 40, 200 10 4, 20 20%, 15% 

The first simulation (termed pessimistic - P), assumes 
that only 210 clusters are distributed in the redshift 
ranges in the following form: < z < 0.5 (100), 
0.5 < z < 1 (70) and 1 < z < 1.5 (40) with ADD statisti- 
cal errors of 15%, 17% and 20%, respectively (see Table 
!)■ 

In the second one (optimistic case - O), 1050 clusters 
were redshift distributed as follows, < z < 0.5 (500), 
0.5 < z < 1 (350) and 1 < z < 1.5 (200) with ADD sta- 
tistical errors of 10%, 12% and 15%, respectively. The 
redshift intervals were partitioned into bins (Az = 0.05) 
with the clusters distributed as shown in Table l[20j. 
Both simulations were carried out by marginalizing over 
the Hq parameter in DA{zi,p) in Eqs. (1) and (5). 
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FIG. 3: Contours on the plane (go,gi) from synthetic Planck 
data as defined in Table 1. In the pessimistic case, 210 P 
(pessimistic projection) and 1050 O (optimistic projection) 
clusters were considered by assuming a random error between 
10% and 20% in the resulting angular diameter distances. 

In Fig. 3 we display the results of our simulations for 
the linear parameterization. The contours correspond 
68%, 95% and 99.7% c.l. for the optimistic (colored in- 
ner contours) and pessimistic case (outer contours), re- 
spectively. Comparing with Fig. 1(a) we see that the 
allowed region was remarkably reduced even in the pes- 
simistic case. This means that angular diameters from 
SZE/X-ray data will become a potent tool for measuring 
cosmological parameters fully independent and competi- 



tive with SNe type la luminosity distances. 

Conclusions. We have shown that the combination of 
Sunyaev-Zeldovich/X-ray data from galaxy clusters is an 
interesting technique for accessing the present accelerat- 
ing stage of the Universe. This result follows from a new 
kinematic approach based on the angular diameter dis- 
tance of galaxy clusters obtained from SZE/X-ray mea- 
surements. By using two different parameterizations, it 
was found that the existence of a transition from a decel- 
erating to an accelerating expansion occurred relatively 
recent (zj ~ 0.3). 

The ability of the ongoing Planck satellite mission to 
constrain the accelerating stage was discussed by sim- 
ulating two realistic samples of angular diameters from 
clusters. The allowed regions in space parameter was sig- 
nificantly constrained for both the pessimistic and opti- 
mist simulations (Figure 3) . The limits on the transition 
redshift derived here reinforces the extreme interest on 
the observational search for obtaining SZE/X-ray data 
from galaxy clusters. 

Finally, we stress that the present results depends nei- 
ther on the validity of general relativity nor the matter- 
energy contents of the Universe and, perhaps, more im- 
portant, they are also independent from SNe type la ob- 
servations. 
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